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Abstract 

In this work we present ab initio calculations of the formation energies 
and stability of different types of multi-vacancies in carbon nanotubes. We 
demonstrate that, as in the case of graphene, the reconstruction of the de- 
fects has drastic effects on the energetics of the tubes. In particular, the 
formation of pentagons eliminates the dangling bonds thus lowering the for- 
mation energy. This competition leads to vacancies having an even number 
of carbon atoms removed to be more stable. Finally the appearance of magic 
numbers indicating more stable defects can be represented by a model for 
the formation energies that is based on the number of dangling bonds of the 
unreconstructed system, the pentagons and the relaxation of the final form 
of the defect formed after the relaxation. 
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1. Introduction 

Nanoscopic systems have attracted significant attention from the scien- 
tific community due to the possibility of designing ever-smaller electronic 
devices. Amongst candidates with greatest potential for application one can 
find carbon-based materials such as carbon nanotubes (CNT)[l|, and more 
recently graphene^]. 
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Intrinsic defects are often seen as the source of deleterious effects in semi- 
conductor materials jij. For instance, vacancies and clusters of vacancies have 
important well known effects on the properties of many semiconductors of 
technological importance such as SzfHJ, GaAsjfjj], SiGep], [s[ and Gej^]. The 
same, however, cannot be clearly stated about carbon based materials. In 
fact, defects in carbon nanotubes can be used as binding sites for different 
types of gaseous species in CNT-based sensors [ill. 1 1 01] - as well as in a new 



family of disordered graphene spintronics devices [12]]. One point is clear: de 



fects can lead to drastic changes in the electronic structure of carbon-based 
systems. Thus if these materials are to be the building blocks of tomorrow's 
electronics, characterizing these defects is of utmost importance. 

In carbon nanotubes one can find many different types of defects ran ging 



from Stone- Wales 13 , to adatoms 16, 15, 14 and vacancies 22, 25, 19, 26 



2Q|, |2J, [271, [18|, [171, |21|, |23j. Recently, Saito et a/.[28j have reported on the 
existence of point defects in graphene, ranging from a single vacancy to an 
octovacancy[29]. The authors also show that the divacancy, the tetravacancy, 
and the hexavacancy, are the most stable defects in a graphene sheet. 

In order to highlight the stability of specific vacancies one introduces the 
concept of magic numbers 28j . These magic numbers indicate the size - the 
number of atoms removed - of a defect that lead to the most stable multi- 
vacancies. Previous studies using positron annihilation in graphite suggest 
that the hexavacancy, V 6 , was the most stable defect [24]] . This stability was 
explained by the dangling bond counting model (DBCM), where the number 
of dangling bonds (DB) of the system, Ndb, decreased as the vacancy be- 
came more stable. This model was also successful in explaining the stability 
of multivacancies in silicon|5j and GaAs|6]. 

In graphene, however, Saito et al. [28], noted that the stability of the 
defect also depends on pentagons formed upon reconstruction of the system. 
Subsequently extending the dangling bond counting model - adding to that 
the effect of the pentagons - the authors proposed the pentagon and dangling 
bond counting model (PDBCM). The PDBCM is based on the average energy 
per DB and per pentagon in all defects considered. The model was then 
used in explaining the stability of vacancies (and consequently the existence 
of magic numbers) in graphene. 

The aim of the present study is twofold. The first one is to determine 
what are the most stable vacancies in carbon nanotubes and whether the 
PDBCM is applicable in the presence of hybridization between ir orbitals 
due to curvature effects. 



2 



Secondly and most importantly we propose a modified model based on 
the PDBCM that uses only two defects for the construction of the model. In 
other words, one does not need to perform calculations for all vacancies to 
determine the parameters for the model; all the parameters - as we will show 
- can be obtained by two preliminary calculations. Most importantly, as it 
will be demonstrated, it contains all the physical ingredients to predict the 
existence of magic numbers in carbon nanotubes, even in the case of large 
curvature effects. 



2. Method 

Ab initio total energy calculations based on density functional theory^' 



31] were performed for different types of vacancies on (5, 5), (7, 7), (9, 9) and 
(10, 10) carbon nanotubes. This smaller-radius nanotubes was chosen in or- 
der to give rise to large curvature effects, as oposed to the plan ar g raphene 
sheet. We used the generalized gradient approximation (GGA)[32[ for the 
exchange and correlation potential within the Perdew-Burke-Ernzerhof ap- 
proach [33j. Our simulations were performed using a plane- wave DFT method 



within the VASP code|35l. I34j| and with ultrasoft pseudopotentials. [36] In all 
our calculations a plane wave energy cutoff of 290 eV and 8 k-points along 
the reciprocal axis of the CNT were used. 

Initially the defects are created by simply removing n (n = 1 . . . 8) carbon 
atoms from a pristine structure containing 8 irreducible unit cells for each 
nanotube[37||. The (5,5) nanotubes with different numbers of carbon atoms 
removed prior to relaxation are shown in figure []Jal — hi) - left hand-side 
panel ( the defects on the other nanotubes are equal). Hereafter we label 
V n the respective vacancy where n carbon atoms have been removed. We 
note that the creation of vacancies leads to the formation of dangling bonds 
- carbon atoms with two-fold coordination instead of the expected three- 
fold one - which are energetically unfavorable. The systems are allowed to 
atomically rearrange using a conjugate gradient method (CG) until the forces 
on all the atoms are lower than 0.02 eV/A. The final relaxed structures are 
shown in figure [T](a2 — hi) - right hand-side panel. One can notice that the 
defects undergo a reconstruction that leads to the formation of pentagons 
and subsequent saturation of the dangling bonds. We note that in the case 
of n even, the number of DBs goes to zero (except for V$). 
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3. Results 



The formation energy, Ef[n] for the V n vacancy is calculated using, 

E f [n] = E r [n] - E p [n] - nfi c , (1) 

where E p is the total energy of the pristine nanotube, E r [n] is the total 
energy of the reconstructed nanotube with a V n defect, n is the number of 
carbon atoms removed from the system and /j,c is the chemical potential of 
1 carbon atom, which is the total energy of the pristine nanotube divided by 
the number of atoms on the sysmtem. 

The total formation energy as a function of vacancy size is shown in figure 
|2}i. From the figure, one might be tempted to say that the Vi is the most 
stable vacancy. However, a more reasonable approach is to compare each 
vacanc y n ormalizing the number of atoms that have been removed from the 



system [38j. In that case we observe from figure© that the hexavacancy has 
the lowest formation energy per C atom removed. We also note that Ef [n] is 
nonmonotonic. Instead, it has local minima for even n. This result is in line 
with previous results for graphene indicating the existence of the so-called 
magic numbers, namely 2, 4 and 6 for carbon nanotubes as well as graphene. 

That, however, is not the full picture. One step further into fully un- 
derstanding the stability of vacancies in CNTs is to determine how they are 
correlated with closely sized vacancies, for instance, whether a V§ will break 
into a single vancancy, V\, and a pentavacany, V5, and so on. 



Hence, in order to address defect stability Saito et a/. [28] proposed two 
quantities associated with distinct dissociation processes. The first one as- 
sumes that a V n — type defect breaks up into a single vacancy, V\, and a 
V n -i — type vacancy, 

V n K-l + V 1 . (2) 

In that manner the first dissociation energy is defined as the energy change 
between initial and final states, 

Di [n] =E f [n-l] + E f [1] - E f [n] . (3) 

For the second case, one would have two V n —type defects that reconstruct 
into a V n _i — type and a V n+ i — type vacancy. In other words, a single vacancy 
breaks away from one of the defects and migrates to another one close by 



4 



2V n ->• K-i + K+i • (4) 
Thus the second dissociation energy considered here is defined as 



D 2 [n] = E f [n + l} + E f [n - 1] - 2E f [n] . (5) 

We finish by highlighting that, following the above definition, the higher 
the value of D\ [n] and D 2 [n], the more stable the defect is. 

In Figure [3](a — b) we observe peaks in the dissociation energy, localized in 



we can con- 



V 2 , V4 and Vq. Following the same idea given by Saito et al. |28 
elude that the vacancies V 2 , V4 and V e are stable in the carbon nanotubes in 
a fashion similar to graphene, corroborating the initial conclusion that there 
are three magic numbers. Furthermore, the curvature for such nanotubes 
is quite large, therefore we can also conclude that these magic numbers are 
valid for graphene and nanotubes in general. 

In order to understand the origin of the stability for these particular 
defects we will focus on the nanotube that has the largest curvature effect 
- the (5,5). First we will look to a model that only takes into account the 
dangling bonds, namely the dangling bond counting model ( DBCM)^2^. 
We therefore need to determine some quantities. The first is the gain in 
energy due to the relaxation of the system, which is defined as the difference 
in energy between the reconstructed system and the energy of the system 
with the unrelaxed defect, E u , 



E re lax [n] = E u [n] - E r [n] . (6) 

The pentagon bond energy is then defined as the ratio between the re- 
laxation energy and the number of pentagons, Np ent , in the relaxed CNT 



Epent M = E re i ax [n] /N pent . (7) 

Finally we also define the energy per dangling bond as the ratio between 
the formation energy of the unrelaxed defect and the number of dangling 
bonds 



5 



E D b [n] = (E u - E p [n] - nji c ) /N DB 



(8) 



In Table [T] we summarize these quantities for each one of the vacancies in 
the present study obtained from our DFT calculations. 
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Table 1: Energies associated with atomic relaxations, E re i ax , pentagons, E pent , 
dangling bonds, Ed B for each type of vacancy. The total number of pentagons, 
N p ent, dangling bonds (prior to, N^B and after, N^f^, reconstruction), and the 
number of sides, N s ides of the reconstructed poligon are also shown. The last 
column indicates the average of some of these quantities over all the defects. 



The first model we analyze only takes into consideration the dangling 
bonds in the system. Thus the formation energy in the so called dangling 
bond counting model is simply given by 

Ef BCM [n] = 2A9N DB (9) 

where N DB is the number of dangling bonds for the unreconstructed system 
and the proportionality factor corresponds to the energy per DB, averaged 
over all possible vacancies as shown in Table [TJ 

The DBCM is a rather simplified model and it does not take into account 
the strong relaxation of the defects. On the other hand the pentagon and 
dangling bond counting model proposed by Saito et al. also considers the 
reconstruction of the vacancies into pentagons. In the PDBCM the formation 
energy is thus obtained by considering also the reconstruction of the vacancies 
into pentagons. Hence the pentagon and dangling bond counting model adds 
a correction to the DBCM to include the gain in energy due to the formation 
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of pentagons. This new term includes the number of pentagons and the 
average pentagon bond energy. Thus the formation energy in the PDBCM 
is written as 

e pdbcm [n] = 2A9N DB - 2.50N Pent (10) 

where Npb is the number of pentagons of the relaxed system. We note that 
the formation energy is now given by a competition between the high-energy 
dangling bonds and the reconstruction of the defects into pentagons to try 
to eliminate as many DBs as possible. 

In figure H] we present the formation energy calculated with both methods, 
DBCM and PDBCM, compared with our DFT results assumed here to be a 
benchmark calculation. 

From figure H] we can extract two pieces of information. First, if we 
consider only the dangling bonds, we obtain a linear relationship over a 
wide range of defects. That is clearly not the profile seen by performing 
the full DFT calculation. It is then clear that the reconstruction plays an 
important role in the formation energy of the defect. Using the PDBCM, 
where not only the dangling bonds are taken into consideration but also the 
number of pentagons formed after relaxation, we can see that the results are 
in reasonable agreement with our DFT calculations. 

We also use both models to fit the dissociation energies and to determine 
whether they reproduce the magic numbers that have been found in our 
DFT calculations. From figure - which shows the dissociation energies as 
a function of defect size - we can conclude again, that the dangling bond 
counting model fails to reproduce the two dissociation energies, whereas the 
model that takes into account the pentagons, reproduces reasonably well 
(specially in the case of D2) our DFT calculations. 

Thus, the PDBCM as proposed by Saito et al, can be used to explain 
the appearance of the magic numbers and the stability of the vacancies. The 
model presented here uses averages over all defects, but even if one uses 
only two defects - the single vacancy and the V2 - would lead to a similar 
result. One downside of the PDBCM, however, is the fact that it cannot 
account for the dips in the formation energy for even-numbered vacancies. By 
construction, the formation energies of these vacancies are always identical to 
the previous odd-numbered ones. Thus, Ef [n], it is not possible to infer the 
existence of the magic numbers. This leads to the conclusion that the model 
is missing an important ingredient. In the light of this problem we propose 
a new model for the formation energy of multivacancies. This model retains 
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the spirit of the PDBCM where DBs and the reconstruction of the defect 
are the main ingredients for the formation energy, but it should also include 
information about the final shape of the defect. In others words, energy is 
gained by forming the pentagons, but in detriment of other bonds that are 
bent after the relaxation is complete. 

Our proposal for the formation energy is based on three main contribu- 
tions. The first includes the contribution due to the dangling bonds while 
the second accounts for the formation of pentagons after relaxation. Both of 
which retain the same spirit of the PDBCM, 

E f [n] = E DB N DB - E P Np + E S (N S -N P ). (11) 

The third and final term accounts for the relaxation energy (E s ) of all bonds 
in the defect which have not been included in the pentagons. 

In order to determine the parameters described above, we choose the 



single vacancy and the trivacancy [39] . Since both show, upon relaxation 



pantagons and dangling bonds. The dangling bond energy is obtained by 
taking the average between two defects. We then use the DFT results for 
V\ and V3 in equation ( FIT]) together with the average dangling bond energy. 
This leads to a system of simultaneous equations if one assumes the DFT 
result to be the correct formation energy in this case. The parameters are 
thus summarized in table [2J 



E DB (eV) E P (eV) E s (eV) 
2.57 -3.28 0.15 



Table 2: Parameters used in our calculations of the formation energies of 
n— vacancies in a (5, 5) carbon nanotube. 

The formation energies for all defects, calculated from equation (1111) are 
depicted in figure |6j One can see that the MPDBCM is, in general, better 
than the PDBCM. In particular it is capable of accounting for the dips on 
the formation energy of even-numbered multivacancies. The only case where 
the MPDBCM performs worse is the tetravacancy. For the V4, the pentagon 
oriented along the tube has much longer bond lenghts when compared to the 
C — C bonds ( 1.62A as oposed to 1.50A). This leads to a discrepancy in the 
formation energy that coincidentally is correct in the PDBCM. 
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The dissociation energies calculated with the MPDBCM (shown in figure 
[7]) are also in very good agreement when compared with our DFT calcula- 
tions. Again they are usually better. 

4. Conclusions 

Hence we have calculated the formation and dissociation energies asso- 
ciated with n-vacancies in carbon nanotubes. We have demostrated that in 
carbon nanotubes with large curvature effects as in graphene the same sta- 
ble multivacancies appear, namely the divacancy, the tetravacancy and the 
hexavacancy. The existance of these magic numbers is corroborated by the 
pentagon and dangling bond counting model which shows that the stability of 
the defects is given by the competition between high energy dangling bonds 
and the reconstruction of the multivacancies. In the cases where the defects 
are most stable we note that the relaxation leads to no dangling bonds left. 

Finally we proposed a modified method based on the PDBCM which also 
takes into account the final shape of the system. It takes into consideration 
both dangling bonds and the reconstruction of the defect and adds to that 
effect of bond streching and bond bending that takes place of reconstruction. 
This way we are able to include one important physical ingredient that was 
missing from previous models. From our calculations one can obtain the 
parameters that fit extremely well the formation and dissociation energies of 
a number of n-vacancies. In particular, our model is capable of accounting 
for magic numbers 2,4 and 6 which are related to the most stable defects. 
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(d1) (d2) 




Figure 1: Unrelaxed (left hand-side panel) and relaxed (right hand-side panel) 
multivacancies: (a) single vacancy, V\\ (b) divacancy V2; (c) trivacancy V3 (d) 
tetravacancy V4; (e) pentavacacy, V5; (f) hexavacancy, Vq; (g) heptavacancy, V7; 
(h) octavacancy, Vs- 
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1 2 3 4 5 6 7 8 

Vacancy Size (n) 

Figure 2: (a) Total formation energy, and (b) formation energy per carbon atom 
removed for a (5, 5), (7, 7), (9, 9) and (10, 10) carbon nanotubes obtained via DFT 
calculations. 




Vacancy Size (n) 

Figure 3: Dissociation energies (a) D\ [n] (equation ([5])), and (b) D2 [n] (equa- 
tion fl5J)) as a function of defect size for a (5,5), (7,7), (9,9) and (10,10) carbon 
nanotubes. 
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Figure 4: Formation energies calculated using the DBCM (blue triangle) and PDBCM 
(red diamonds) . The DFT calculations using equation Q] (black square) are used as a 
benchmark for the quality of both models 
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Figure 5: Dissociation energies a) D\ [n] (equation ©), and b) D 2 [n] (equation ([5])) for 
the dangling bond couting model and the pentagon and dangling bond counting model. 
We also present our DFT results for comparsion. 
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Figure 6: Formation energies calculated with the Modified-PDBCM (equation 
(jlip ). and with the PDBCM (equation (jlOp ). Our benchmark DFT calculations 
(black square) are also shown for the sake of comparison. 
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Figure 7: Dissociation energies as a function of the number of vacancies, a) D\ [n] and b) 

D 2 [n] 
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